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Abstract 



Collisionless plasmas are ubiquitous in the Universe. In general full kinetic equations are 
needed for a correct description of the physics in such environments. Fluid approximations 
are, at some special cases, valid. In this work we present the basic equations under such 
an approximation and the implementation of them in a magnetohydrodynamical code. The 
numerical simulations are applied for intracluster medium physical properties, i.e. low density 
but high temperature plasma. We show that, under specific conditions firehose and mirror 
instabilities may be present. The side effects are the conversion of thermal to magnetic/kinetic 
energies and the generation of small scale structures (~ lkpc). We study the changes in 
the statistics of Faraday rotation measurement, and of the synchrotron polarization, when 
compared to a standard MHD simulation for the ICM. We conclude that KMHD numerical 
simulations are mandatory in order to fully understand the statistics of the ICM emission. 



1 Introduction 



The physical properties of the interstellar medium, such as mean magnetic field (B ~ /iG), 
temperature (T ~ 20 — 10 4 K) and density (10 — 10 4 cm -3 ), result in mean-free paths (A m f p ~ 
10AU) much smaller than typical dimensions of the system (L ~ lOOpc). It means that 
the ISM plasma may be well described under fluid approximation. Also, since the collision 
frequency of ions is typically much larger than their cyclotron frequency (t>n 3> Sli), we may 
^ ■ also assume the momentum distribution to be isotropic (Maxwellian) . 

In the intracluster medium, on the other hand, with (B ~ 0.1/iG), temperature (T ~ 
10 8 K) and density (10 -2 cm -3 ), the collision frequency of ions is typically much smaller than 
their cyclotron frequency [y^ <C though its mean-free path (A m f p ~ lkpc) is still smaller 
than the typical dimensions of the system (L ~ 0.1 — IMpc). Therefore, even though the ICM 
being properly described under fluid approximation, its momentum distribution will no longer 
be isotropic. These plasmas are often called gyrotropic and the momentum distribution of 
particles is usually described by a bi-Maxwellian function. One for the motions parallel to 
the magnetic field, and another for the perpendicular. In other words, a gyrotropic plasma 
presents two temperatures (Tj| and Tj_), with respect to the local magnetic field. 

In this case, some dynamical properties of the plasma may change. The evolution of the 
turbulent cascade will be different, when compared to that of a standard MHD turbulence, 
for instance. Also, other kinetic effects, such as instabilities may take place. If these effects 
are dominant, the understanding of the physical properties of the ICM may be wrong. This 
issue is briefly discussed in this work. 
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2 The KMHD model 

In the fluid approximation, a gyrotropic plasma can be described by the CGL-magnetohydrodynamic 
(Chew et al. 1956) equations expressed as follows (see also Kowal, Falceta-Gongalves & Lazar- 
ian 2011): 
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where p and v are the plasma density and velocity, respectively, B is the magnetic field, 
P = p±I + (p\\ —p±)bb is the pressure tensor, b = B/|B| is the unit vector along the magnetic 
field, p|| and p± are the pressure components parallel and perpendicular to b, respectively, 
and f represents the turbulent forcing term. Further details on the code may be obtained in 
Falceta-Gongalves et al. (2010a, b), Falceta-Gongalves, Lazarian & Houde (2010). 

The above set of equations is closed by the determination of the parallel and perpendicular 
pressures. Here, we make use of the double-polytropic equations, in conservative form, as: 



a(P f 7l) + V.M^ v) = o, (4) 



— g t J - + V-(p ll (B/p)^- 1 v) = 0, 



(5) 



where 7^ and 7|| are the polytropic exponents for the perpendicular and parallel pressures, 
respectively. 

In this work we focus on the double-isothermal case, i.e. the pressure ratio cj_/c|| is 
constant (Hasegawa 1969). In this way we can observe the role of the instabilities during the 
whole simulation. Typically, the firehose and mirror instabilities grow fast, at small scales 
mainly (see Figure 1), and the pressure anisotropy tends to be reduced very quickly. Several 
natural processes, such as scattering of cosmic rays, or anisotropic thermal conduction, may 
lead to the unstable regimes, but we avoided dealing with these complications and focus only 
on the unstable regime. 

In Figure 1 we show the turbulent velocity spectrum of the simulated plasmas at different 
unstable regimes, firehose and mirror. Compared to the standard MHD case, it is clear the 
increase in the total power at small scales due to the instabilities. The mirror unstable case 
is more prominent since it acts at the magnetosonic modes (compressible), while firehose at 
the Alfvenic mode. 



3 The KMHD effects in the ICM 



Among the main observational methods to infer the ICM magnetic fields are the synchrotron 
polarization and the Faraday rotation measurements (e.g. Giovannini, Tordi & Ferreti 1999, 
Ensslin & Vogt 2003, Govoni & Feretti 2004, Guidetti et al. 2010 and others). 
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Figure 1: Velocity power spectrum obtained for the KMHD firehose- unstable model (left), 
with c_i_/c|| = 0.5, and for the KMHD mirror-unstable case (right) with c±/c\\ = 2.0. MHD 
case is shown as dashed line for comparison (extracted from Kowal, Falceta-Gongalves & 
Lazarian 2011). 

The intrinsic large polarization degree of the synchrotron emission makes it very useful 
in studying the ICM magnetic fields. Polarization maps reveal the field line structure in 
the plane of sky. Complex, non-uniform, distributions reveal a kinetic dominated plasma, 
where -E> 2 ky < pSv 2 ~ SB 2 . Very uniform polarization vectors, on the other hand, reveal 
strong magnetic fields with B 2 ky >> pSv 2 (Falceta-Gongalves, Lazarian and Kowal 2008). 
Unfortunately, this technique does not tell much about the magnetic field component parallel 
to the line of sight. In this particular case, the Faraday rotation effect may be useful. 

Most of the radio sources in clusters of galaxies are embedded in magnetized, hot and 
low density plasmas that account for the frequency dependent refraction index, which is 
responsible for the Faraday rotation. The polarization angle of the observed radiation depends 
on the frequency as follows: 

0obs(A) = <£o + A 2 x RM, where (6) 

pL(kpc) 

RM ~ 812 / n e {cmr s )B\\{iiG)dl (rad nT 2 ). (7) 

Jo 

Once the distribution of matter (n e ) is known - from X-rays maps, for example - it is 
possible to extract information about Bu . Combined to the polarization maps, the rotation 
measure (RM) is a good tracer of the magnetic field in the ICM, and its coupling with the 
local plasma. 

Galactic motions within the cluster medium is expected to occur at scales of ~ IMpc, 
with dispersion of velocities of order of 1000 km s _1 . This kinetic energy density is about 
an order of magnitude larger than the energy density of the average uniform magnetic field 
component. Therefore, the turbulence in the ICM tends to be responsible for changes in the 
structure of the field lines as well. 

The statistical study of the RM observed in clusters of galaxies has also been used in 
order to understand its structure and coupling with the turbulent medium (Ensslin & Vogt 
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Figure 2: Up: Rotation measure calculated for the MHD model (left), with Bq = 0.1/iG, L = 
IMpc and ri = 0.01cm" 3 , and for the KMHD mirror-unstable case (right) with ci/cy = 2.0. 
Bottom: power spectrum of the spatial distribution of the rotation measure calculated for 
the MHD model (left), with B„ = 0.1/iG, L = IMpc and n„ = 0.01cm- 3 , and for the KMHD 
mirror-unstable case (right) with c±/cu = 2.0. 

2003). Usually, the observed RM maps are compared to cosmological simulations that are, in 
essence, MHD. In a more realistic case, under KMHD regime, the structure of the magnetic 
field and its role on the propagation of turbulent modes may change. This may be obtained 
from the observed RM maps. In Figure 2 we show the polarization maps obtained for the two 
regimes (left:MHD, right:KMHD - mirror unstable), considering L — lOOkpc, Bo = 0.1/iG 
and < n e >= 0.005. The difference in structure is clear. The KMHD mirror-unstable case 
shows a power increase at small scales, as seen in the power spectra (Figure 2, bottom panels). 
Also, the rotation measure in Figure 2 is larger for the KMHD case. That occurs because 
the mean magnetic field is in the plane of sky. The stronger perturbations in the KMHD 
case bend the field lines towards the LOS direction, increasing the RM. This difference will 
be reduced as the angle between the mean field and the plane of sky increases. 

The polarization maps for both cases are shown in Figure 3. The maximum polarization 
degree obtained for the MHD case is 55%. The intrinsic synchrotron polarization, around 
70%, is reduced due to the turbulent structure of the ICM. The increase of the random 
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Figure 3: Polarization maps of the synchrotron emission calculated for the MHD model (left), 
with Bq = 0.1/iG, L = IMpc and uq = 0.01cm -3 , and for the KMHD firehose-unstable case 
(right) with c±/c\\ = 0.5. 



component of the magnetic field increases the dispersion of the polarization angle (0) and 
reduces the polarization degree. For the KMHD case we see two different situations. For the 
mirror-unstable regime the instabilities produce compressions that are not in the Alfvenic 
modes, so they do not increase the random component of the field. In the firehose-unstable 
regime, on contrary, the instability operates at the Alfvenic modes and result in an increase in 
the random component of the field. As a result, as shown in Figure 3 (right), the polarization 
degree decreases to ~ 33%. 

4 Conclusions 

Main results are: 

• kinetic effects may change the dynamics of the ICM plasma; 

• depending on the regime, fast increase of magnetic field may take place; 

• as well as acceleration (fluid motions), both at small scales (~ lkpc); 

• the KMHD instabilities result in modifications of statistics of RM and synchrotron polar- 
ization. 

In the double adiabatic case, the differences between the standard MHD approach and 
the one presented in this work will, obviously, be less prominent since one of the effects 
of the instabilities is to decrease pressure anisotropy. However, we believe that most of 
the conclusions, mostly regarding fast B-growth and modifications in the statistics of RM 
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and synchrotron polarization remain. This will be tested in the near future with proper 
simulations. 
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